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Erratum 

"Energy  Principle  V/lth  Specific  Inductance" 

By  Glenn  Bateman 

Page  16  bottom:   Delete  "In  particular,  we  see  from  the  ideal 

MHD  equation  of  motion..."  to  end  of 
paragraph. 

Replace  with:   "For  this  minimization  we  must  demand  that 
each  element  of  the  plasma  is  able  to  communicate  with  the 
rest  of  the  flux  surface.   In  particular,  we  exclude  from 
this  theory  those  systems  in  which  all  the  field  lines  are 
closed. " 

Page  17: the  equation  after  Eq.  (28):  Delete  "2ir". 

Page  20  top:   Delete  "The  variation  of  flux  and  inductance 

separately. . . "  to  the  end  of  paragraph. 

Page  23  top:   Delete  "However,  as  the  shells  change  shape,...' 

to  end  of  paragraph.   Replace  with:   "It  can 
be  shown  that  individual  fluid  elements  also  do  not  compress 
or  expand  during  their  motion  within  a  shell  as  they  reach 
equilibrium. " 

Page  35  end  of  first  sentence:   Delete:  "...before  comparing 

states  with  a  variational  analysis." 
Replace  with:   "...in  successive  stages  within  the  class  of 
accessible  states." 

the  specific  Inductance  matrix.   The  latter  depends  upon  only 
the  geometry  of  the  surfaces.   Stability  conditions  are 
derived  by  varying  the  flux  profile  and  the  geometry  (specific 
Inductance).   Also,  a  nonlinear  condition  is  derived  for  the 
special  case  of  finite  geometric  distortions  of  the  flux  sur- 
faces holding  the  flux  profiles  fixed.   Exact  results  for 
the  specific  Inductance  are  calculated  for  straight  cylinders, 
axlsymmetric  toruses,  and  helically  symmetric  configurations, 
all  with  arbitrary  cross  sections. 
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Abstract 

The  determination  of  stability  for  a  scalar  pressure, 
toroidal  plasma  has  been  formulated  by  means  of  an  energy 
variational  principle  in  which  the  geometric  problem  of 
variation  of  the  flux  surface  shapes  is  separated  from  the 
variation  of  the  flux  profile.   By  means  of  successive  mini- 
mizations, holding  the  surface  shapes  and  the  flux  profiles 
fixed,  the  potential  energy  is  written  as  a  one-dimensional 
integral  over  a  function  of  pressure,  differential  flux,  and 
the  specific  inductance  matrix.   The  latter  depends  upon  only 
the  geometry  of  the  surfaces.   Stability  conditions  are 
derived  by  varying  the  flux  profile  and  the  geometry  (specific 
inductance).   Also,  a  nonlinear  condition  is  derived  for  the 
special  case  of  finite  geometric  distortions  of  the  flux  sur- 
faces holding  the  flux  profiles  fixed.   Exact  results  for 
the  specific  Inductance  are  calculated  for  straight  cylinders, 
axisymmetric  toruses,  and  helically  symmetric  configurations, 
all  with  arbitrary  cross  sections. 


-Ill- 


I.    Introduction 

The  energy  principle  for  a  scalar  pressure  plasma  can  be 
reduced,  after  aujxillary  minimizations,  to  a  form  involving 
the  geometry  only  through  a  specific  inductance  matrix. 

The  motivation  for  using  the  magnetic  inductance  is 
illustrated  by  the  following  model.   Consider  the  vacuum  magnetic 

field  enclosed  between  two  nested  toroidal  conductors.   Blank, 

2 
Friedrichs  and  Grad  show  that  the  magnetic  energy  can  be  written 

as  the  product  of  the  total  currents  and  the  total  magnetic 
fluxes  (Hodges'  theorem).   Furthermore,  the  fluxes  can  be 
written  as  a  linear  combination  of  the  total  currents.    '    The 
factors  of  proportionality,  called  the  inductance  matrix,  depend 
upon  only  the  geometry  of  the  conductors  and  their  relative 
position.   Therefore,  the  behavior  of  the  energy  can  be  investi- 
gated by  separately  varying  the  geometry  (Inductance)  or  the 
total  currents.   Once  the  inductance  is  known,  the  energy  can  be 
calculated  from  the  total  currents  without  first  calculating  the 
distribution  of  current  density  on  the  conductors  or  computing 
the  magnetic  field. 

In  this  paper,  the  above  ideas  are  modified  to  apply  to 
toroidal  systems  (plasmas)  with  distributed  current.   Given  a 
set  of  simply  nested  toroidal  flux  surfaces  and  a  pair  of  total 
azimuthal  and  longitudinal  currents  for  each  surface,  there 

exists  one  magnetic  field  configuration  which  minimizes  the 

5  6 
magnetic  energy  on  each  surface.^'    For  this  magnetic  field,  we 

show  that  the  total  currents  can  be  written  as  a  linear  combination 
of  the  fluxes  trapped  between  differentially  close  flux  surfaces. 
The  four  coefficients  in  the  linear  combination  are  called  the 
specific  inductance  matrix.   For  each  flux  surface, the  specific 
Inductance  depends  upon  only  the  geometry  of  the  surface  and  its 
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geometric  relation  to  differentially  close  neighboring  surfaces. 
This  is  proved  in  section  II-G.   Grad  and  Rubirr  and  Kruskal  and 
KulsrucT  show  that  this  minimized  magnetic  energy  reduces  to  a 
one  dimensional  integral  over  the  product  of  total  currents  and 
differential  fluxes.   Using  the  specific  Inductance,  the  integrand 
is  factored  into  geometric  and  profile  dependent  matricies. 
We  may  then  vary  the  energy  by  varying  the  geometry,  holding  the 
flux  profile  fixed  as  a  function  of  volume,  or  by  varying  the 
flux  profile  while  holding  the  geometry  of  the  surfaces  fixed. 
If  we  use  an  ideal  magnetohydrodynamic  model  and  minimize 
the  magnetic  and  thermodynamic  energy  on  each  surface,  the  entire 
potential  energy  reduces  to  a  one  dimensional  integration.   If 
the  appropriate  constraints  are  used  to  assure  that  states  are 
accessible  to  each  other  by  a  possible  motion  of  the  plasma,  it 
can  be  shoi\/rr'   that  the  state  of  extremum  potential  energy 
corresponds  to  the  classical  definition  of  static  equilibrian. 
We  call  the  system  unstable  if  the  potential  energy,  with  the 
magnetic  and  thermodynamic  energy  minimized  on  each  surface,  is 
lower  in  any  other  accessible  state  than  in  the  equilibrium  being 

considered.   In  section  III-E  we  show  that  this  definition 

7 
corresponds  to  the  definition  used  in  the  classical  energy  principle.' 

In  order  to  be  able  to  use  the  variational  formulations  of 
stability  derived  in  sections  III-B,  C,  and  D,  we  must  be  able  to 
calculate  the  specific  inductance  for  both  the  equilibrium  and  all 
the  perturbed  states  to  be  considered.   In  s'ection  II-C,  a  number 
of  formal  representations  for  the  specific  inductance  are  derived. 
In  section  IV,  explicit  results  are  found  (in  terms  of  one  integra- 
tion) for  the  specific  inductance  of  straight  cylindrical,  axi- 
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symmetric  toroidal,  and  helically  symmetric  surfaces,  all  with 
arbitrary  cross  sections.   The  Integration  can  be  performed 
analytically  in  the  first  two  cases  when  the  cross  sections 
are  circular  with  arbitrarily  displaced  centers. 

II.   Magnetic  Energy  in  Terms  of  Specific  Inductance 
A.   Flux  Surfaces 

In  this  section  we  shall  characterize  the  class  of  magnetic 
fields  which  are  compatible  with  a  given  set  of  flux  surfaces. 
Setting  aside  those  problems  associated  with  closed  field  lines  (such 
as  island  structure),  we  shall  assume  that  these  surfaces  exist 
and  are  nested  in  a  single  set  around  one  magnetic  axis.   Each 
surface  is  labeled  by  its  volume  V  and  is  represented  by  a  scalar 
function  V  =  V{0      9      r) .      When  an  explicit  representation  is 
needed,  we  use  a  right  handed  coordinate  system  (0  ,0,-,,r)  with 
the  same  topology  as  the  flux  surfaces  as  illustrated  in  Fig.  1; 
9      and  0p  increase  by  unity  the  long  and  the  short  way  around  the 
torus  respectively. 

On  each  toroidal  surface  there  are  two  homology  classes  of 
closed  line  contours  and  cut  surfaces.   We  use  C-,  and  Cp  to 
represent  closed  contours  the  long  and  the  short  way  around 
respectively  and  S-,  and  Sp  to  represent  a  reciprocal  basis  of 
cut  surfaces  which  are  bounded  by  contours  ~Cp  and  +C,  respectively 
(note  the  exchange  of  subscripts).   These  conventions  are  illustrated 
in  Fig.  1  for  toroidal  surfaces  but  they  may  also  be  applied  to 
cylindrical  surfaces  with  periodic  end  conditions.   Any  closed 
contour  on  the  torus  is  homologous  to  integer  multiples  of  the 
representative  contours  C,  and  Cp.   Blank,  Friedrichs  and  Grad 
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discuss  these  concepts  in  more  generality. 

The  magnetic  flux  through  any  cut  surface  S.  is  defined 

by 

^1  ^  r   ds.B  .  (1) 

^^i 

From  Fig.  1  we  see  that  V'-.  is  the  longitudinal  flux  within  the 
flux  surface  while  'A^  i^  "^^e  azimuthal  flux  outside  the  surface. 
^  fl^x  surface  S,,  is  characterized  by  the  property  that  each  flux 
is  independent  of  the  choice  among  homologous  cut  surfaces.   Using 
only  the  divergence-free  property  of  magnetic  fields 


V-B  =  0  (2) 

it  can  be  shown  that  those  surfaces  which  are  ergodically  covered 
by  non-closed  field  lines  are  in  fact  flux  surfaces.   The  total 
current  through  any  closed  contour  C.  is  defined  by  Ampere's  law 


^  ^J^     d|.B  .  (3) 


"i 

Here  I-,  is  the  total  azimuthal  current  outside  and  Ip  is  the  total 
longitudinal  current  within  a  current  surface.   A  current  surface 
has  the  property  that  each  total  current  is  independent  of  the 
choice  among  homologous  closed  contours  on  that  surface. 

Any  magnetic  field  with  a  given  set  of  flux  surfaces  may  be 
characterized  by  a  scalar  potential  in  the  following  way.   From 
Eq.  (2)  and  the  divergence  theorem  it  can  be  shown  that  the 
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magnetic  field  lines  must  everywhere  lie  in  the  flux  surfaces 

B-VV  =  0  .  (4) 

Equation  (4)  implies  that  there  exists  a  vector  field  w  such 
that 

B  =  VV  X  w  . 

When  this  characterization  of  the  magnetic  field  is  inserted 
into  Equation  (2)  we  find 

VV-V  X  w  =  0  .  (5) 

Without  loss  of  generality  i^/e  may  omit  the  normal  component  of  w 

w-VV  =  0  (6) 

It  can  be  shown  that  any  vector  field  satisfying  Equations  (5) 
and  (6)  may  he  expressed  as  a  surface  gradient  of  a  potential 
function 


w  =  V(2)v  =-Vv-VV^l^  .  (7) 


(A  surface  gradient  is  a  gradient  with  the  normal  component 
subtracted  off.)  Hence,  with  only  the  assumption  that  flux 
surfaces  exist,  we  may  write  any  divergence-free  field  B  in 
terms  of  a  scalar  potential  v 
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B  =  VV  X  v(^)v  .  (8) 


Such  characterizations  of  the  magnetic  field  for  equilibrium 
plasmas  are  discussed  by  Grad  and  Rubin  ^  and  by  Kruskal  and 
Kulsrud  .   Although  we  are  not  discussing  equilibria  we  shall  often 
refer  to  their  calculations  when  the  details  follow  in  the  same 
way.   The  potential  v  is  not  single  valued  on  the  torus;  in  both 
of  the  above  references,  the  period  of  v  is  calculated  to  be 


(f    di'V^^U  =  -ip   (V) 


^1 

(9) 


^2 


where 


^  (V)  =  d^  /dV  . 
J        J 


These  periods  may  be  used  to  write  an  explicit  representation 
for  the  potential  and  for  the  magnetic  field 

(10) 
B  =  i^  VV  X   v^^^e^  -  f^   vv  X  v^^^e^  +^v  x  V^^^A  . 

Note  that  A  is  a  single  valued  function  on  the  torus  and  that  it 
depends  upon  the  choice  of  coordinates  (we  shall  see  later  that  it 
can  be  made  constant  on  each  surface  for  the  appropriate  choice 
of  coordinates) .   For  fields  associated  with  a  given  pair  of  fluxes, 
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A  characterizes  the  distribution  of  the  field  over  each  flux 
surface. 

B.   Magnetic  Energy  (Hodges  '  Theorem) 

Up  to  now  we  have  considered  the  class  of  all  possible 
magnetic  fields  associated  with  a  given  set  of  flux  surfaces  and 
with  a  given  pair  of  fluxes  {ip  ,'ip^)    for  each  surface.   We  now 

consider  only  that  magnetic  field  configuration  from  the  above 

r  3   12 
class  which  minimizes  the  magnetic  energy  /  d.  x  — b  .   If  we  use 

the  representation  of  the  magnetic  field  given  by  equation  (10) 

we  see  that  the  single  -  valued  function  A  is  still  free  when  the 

functions  V(x),7//  (v)  and  ^p(^)  ^^^   specified  and  when  any  specific 

choice  is  made  for  the  coordinates  9-,  and  0p.   If  we  extremize  the 

magnetic  energy  by  varying  A  (see  Eqs .  (H6)  and  (H7)  in  Kruskal 

and  Kulsrud  or  Equation  (25)  in  Grad  and  Rubin  ^  we  find 

VV.V  X  B  =  0  .  (11) 

As  described  in  the  last  section,  any  vector  field  which  satisfies 
equations  of  the  form  (11)  and  (4)  may  be  written  as  the  surface 
gradient  [Eq.  (7) ]  of  a  potential  function 

B  =  V^^)*  .  (12) 

The  potential  <J>  is  not  single  valued  on  the  torus;  according  to 
Ampere's  law  [Eq.  (2)],  the  period  of  *  is  the  total  current 
passing  through  the  contour  C. 
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/  d^.V^^U  =  I.  (13) 

Jc.    ~  ^ 


c 

1 

Since  contour  integrals  over  a  surface  gradient  are  the  same 
for  any  choice  of  homologous  contours,  Eq.  (13)  implies  that 
current  and  fliix  surfaces  coincide  and  that  the  total  currents 
are  surface  quantities  I.  =  I.(V). 

From  the  second  variation  of  magnetic  energy  with  respect 
to  A  it  follows  that  the  extremum  is  always  a  local  minimum  of 
the  magnetic  energy  on  each  surface.   All  higher  order  variations 
are  zero.   This  suggests  that  the  extremum  is  an  absolute  minimum. 
This  property  can  be  shown  explicitly  by  using  the  special  coordinate 
system  to  be  defined  by  Eqs .  (21)  and  (23).   It  will  be  shown  that 
the  minimum  is  attained  when  A  is  a  surface  quantity  and  therefore 
drops  out  of  the  representation  for  B. 

The  magnetic  energy,  so  minimized,  may  now  be  written  completely 

in  terms  of  the  total  currents  and  the  differential  fluxes.   This  ap- 
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plication  of  Hodges '  theorem  may  be  demonstrated  explicitly  '   by 

using  Eq.  (10)  and  (11)  to  write 


B^   =  ^V'^xlie^ij/^-  e^Tp^+\]B]     . 


If  the  flux  surfaces  are  labeled  so  that  V  ==  0  is  on  the  magnetic 
axis  and  V  =  C  is  the  volume  of  a  rigid,  perfectly  conducting, 
toroidal  outer  boundary,  the  integration  over  volume  may  be 
written 

/d^x...  =  r  dV  r  ds  ^~r  ....  (l4) 
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We  then  use  Stokes'  theorem  and  the  fact  that  VV"/|VV|  is  the 
unit  normal  to  each  flux  surface  to  write 


M    ,         ^       

mln     '~  "^  1=1   '-^  o 


=  r   d?c   |b2  =  1  XZ    r  dV  I.(V)^.(V)  (15) 


Note  that  the  single- valued  function  A  does  not  contribute;  the 
minimum  magnetic  energy  M  depends  upon  only  surface  quantities. 
Eq.  (15)  is  true  only  when  flux  and  current  surfaces  colncide-- 
only  when  no  current  (  j  s  V  x  B)  passes  through  flux  surfaces. 

C.   Specific  Inductance 

For  any  system  with  simply  nested  flux  surfaces  and 
coincident  current  surfaces  we  now  assert  that  there  is  a  linear 
relation  between  the  total  currents  and  the  differential  fluxes 

2 

ii(v)  =ZI  Ai.(v)^.(v)  (16) 

j=i   ^        ^ 

where  A. .(V)  depends  upon  only  the  geometry  of  the  flux  surfaces 

-'-J 

(is  independent  of  1.  and  v^.).   We  shall  refer  toA.  .  as  the 
specific  inductance  matrix.   If  we  insert  this  linear  relation 
into  Eq.  (15)  we  derive  an  expression  for  the  minimized  magnetic 
energy  in  which  the  geometry  and  flux  profiles  are  separately 
factored 

r   d^x^B^  =\^£ZIi   f    dV  f    {V)Ii^    {V)^P JV)     .        (17) 
'-'min  i,  j=l  ^  o  "^    ^ 

This  formulation  is  due  to  Harold  Grad. 

For  vacuum  fields,  the  specific  inductance  matrix  is  the 
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reciprocal  of  the  derivative  with  respect  to  volume  of  the 
ordinary  inductance  matrix.   (A-  •  should  rigorously  be  called 
the  specific  susceptance  matrix  and  its  reciprocal  called  the 
inductance.   Only  the  more  well-known  term  inductance  will  be 
used  here.   No  confusion  arises  because  there  is  no  need  to 
Invert  the  matrix.)   The  specific  inductance,  relating  the 
total  currents  to  the  fluxes  between  differentially  close  flux 
surfaces,  can  be  applied  to  systems  with  distributed  current 
densities  (provided  the  current  density  lies  in  the  flux  sur- 
faces), while  the  ordinary  inductance,  relating  total  currents 
to  total  fluxes,  is  properly  applied  only  to  vacuum  magnetic 
fields. 

We  shall  derive  three  representations  for  A-  -(V)  and  prove 
that  it  is  a  positive  definite,  symmetric  matrix  which  depends 
upon  only  the  geometry  of  the  flux  surfaces. 

The  most  elegant  derivation  of  the  specific  inductance,  due  to 
Grad  ,  relies  on  a  uniqueness  theorem  for  the  potentials.   Accord- 
ing to  Eqs .  (8)  and  (12),  the  magnetic  field  is  characterized  by 
two  potential  functions  v  and  *>  which  are  the  solutions  to  the 
following  elliptic  system  of  equations 

VV  X  v(^)v  =  v(^U  .  (18) 

By  the  nature  of  the  surface  gradient  [Eq.  (7) ],  Eq.  (18) 
represents  only  two  equations  (one  for  each  component  lying 
in  the  flux  surface)  for  the  two  unknowns  v  and  <t>.   ¥e  use  the 
lemma  '^   that  the  solutions  of  Eq.  (18)  are  uniquely  determined 
on  a  flux  surface  given  only  the  two  periods  of  ^,    which  are  the 
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total  currents  I..   Now  Eq.  (l8)  may  be  rearranged,  "by  crossing 
It  with  VV 


|vv| 

By  the  same  lemma,  the  solutions  to  Eq.  (19)  are  uniquely  determined 
on  a  flux  surface  (provided  |vv|  ^   0)  given  only  the  periods  of  v, 
which  are  the  differential  fluxes  d^ ./dV.   Since  the  same  solutions 
to  this  set  of  linear  equations  are  uniquely  determined  given 
either  set  of  periods,  the  periods  must  "be  linearly  dependent,  as 
in  Eq.  (l6).   The  coefficients  A. .  in  this  linear  relation  be- 
tween  the  periods  must  depend  upon  only  the  common  given  infor- 
mation, which  is  the  geometry  of  the  flux  surfaces. 

In  order  to  compute  the  specific  inductance  matrix  A. .  we 
may  use  any  solutions  v  and  0  of  Eq.  (l8),  compute  their  periods, 
and  insert  into  Eq.  (l6).   If  the  off-diagonal  terms  are  present 
in  A . .,  two  sets  of  solutions  for  v  and  0  with  distinct  periods 
are  needed  in  order  to  invert  the  matrix  Eq.  (l6). 

An  equivalent  but  more  intuitively  appealing  derivation 
follows  if  we  insert  the  characterization  of  the  magnetic  field 
given  by  Eq.  (10)  into  the  condition  that  the  current  density 
lies  in  the  flux  surfaces  [Eq.  (11)].   This  gives  an  elliptic 
equation  for  the  single  valued  function  a(9-,,0q,V") 


.  .  VV  X  VS.- VV  X  V0  .  Afl  n 

+  XI.  ^V.V  X  (VV  X  V0  )[(-l)J+l^    (V)  +>,   ]  =  0 
J  J  J  +x        *J  . 
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(20) 


Suppose  we  choose  a  coordinate  system  such  that 


VV.V  X  (VV  X  V0.)  =0  (21) 


We  are  then  left  with  only  the  second  order  terms  in  the  equation 
for  A 

H.j  Wx  V0..VVX  Ve.A,_Q_  =0  . 

For  the  solution  to  be  a  single-valued  function  on  each  toroidal 
flux  surface  it  must  be  doubly  periodic  in  (99).      By  the  max- 
imal principle,    the  solution  cannot  assume  either  a  local  max- 
imum or  minimum  within  any  boundary  contours  prescribed  on  the 
surface.   Since  any  closed  contours  on  the  torus  may  be  used 
for  the  boundary,  it  follows  that  A  must  be  constant  over  the 
entire  flux  surface  and  therefore  does  not  contribute  to  B  in 
Eq.  (10).   Now  if  we  use  Eq.  (10)  and  Eq.  (3)  to  compute  the 
total  current  I,  we  find 


Aj_j(V)  =  (-l)J  ^   d^-VVx  Ve  ^^  (22) 

^i 

When  using  Eq.  (22)  to  compute  the  specific  inductance  matrix, 
the  coordinates  0.^  and  0p  must  satisfy  Eq.  (21),  a  purely 
geometric  condition,  and  they  must  have  unit  periods 


d£'V9  .   =  5.  .  .  (2^) 

C.  ~   J    ^J 


-12- 


This  is  a  constructive  derivation  which  explicitly  shows  that 
the  specific  inductance  is  a  purely  geometric  entity. 

Note  that  Eq.  (21)  implies  that  the  vector  field  VV  x  V0 . 

J 

must  be  a  surface  gradient  of  some  potential 

VV  X  v^^^e .  =  V^^U'.  . 
3  3 

Hence,  the  coordinates  0      and  0p  may  be  thought  of  as  the  particular 
choices  of  the  magnetic  potential  v  with  normalizations  given  by 
Eq.  (23).   The  lemma,  referred  to  before,  that  solutions  to 
equations  of  this  kind  are  determined  by  their  periods  uniquely 
up  to  a  constant  on  each  surface  can  be  proved  by  the  methods  used 
here.   Suppose  distinct  solutions  with  the  same  periods  existed. 
They  would  differ  by  a  single  valued  function  which  would  satisfy 
the  homogeneous  part  of  Eq.  (20).   The  maximal  principle    is 
then  be  used  to  show  that  this  function  is  constant  on  each  surface. 

o 

It  should  be  pointed  out   that  the  coordinate  system  just 
derived  {{Y,  9   ,0^)    satisfying  Eqs .  (21)  and  (23))  provides  an 

Q 

alternative  to  Hamada  coordinates.    Hamada  coordinates  exist  only 
in  equilibrium  --  when  j  x  B  =  Vp  for  j  =  V  x  B  and  p(V)  a  surface 
quantity.   Our  coordinate  system  exists  in  or  out  of  equilibrium 
provided  the  magnetic  energy  is  minimized  within  each  surface.   Even 
in  equilibrium,  our  coordinate  system  is  equivalent  to  Hamada 

o 

coordinates  if  and  only  if 

VV'V  X  j  -  0  .  (24) 


In  general,  Hamada  coordinates  are  functionals  of  ip •    (i  =  1  or  2) 


-13- 


while  our  coordinate  system  depends  upon  only  the  shape  of  the 
flux  surfaces.   The  field  lines  of  both  B  and  j  are  straight  in 
Hamada  coordinates.   Our  coordinate  system  is  one  of  the  possible 
coordinates  in  which  only  the  B  lines  are  straight. 

If  the  metric  elements  g.  .  =  e..e.,  e.  =  (-l)'^VV  x  ye^'^^/B 
D  =  V0^  X  V9  'VV,  are  known  for  the  coordinates  satisfying  Eqs .  (21) 


and  (23),  it  is  useful  to  use  a  representation  for  /[.  .   derived  by 

8  ^^ 

Gunther  Spies.    Covariant  and  contravariant  components  are  used 

to  express  the  square  of  the  B  field  in  the  magnetic  energy  Eq.  (17) 


B   =  B.B"^.   B-^  =  ^-^D,   B.  =  g.  .B^  . 
1   »  >   1   ^ij 


After  using  Eq.  (1^),  we  may  identify  the  specific  inductance 
factor  in  Eq.  (17)  with  the  expression 

^IJ  =  i   W  Slj"'  '  //^«'«'  8,jD  .  (25) 

From  the  well  known  properties  of  the  metric  tensor,  it  follows 

that  A. .  is  symmetric  and  positive  definite.   Note  that  the 
-^  J 

equation  for  the  coordinates  [Eq.  (21) ]  implies 

(gljD),2  =  (g2.D),l 

so  that  only  one  integration  is  needed  in  Eq.  (25). 
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III.  Formal  Variation  of  Potential  Energy 
A.  Constraints 

We  shall  now  use  a  simple,  scalar  pressure,  magnetohydrodynamic 
model  to  investigate  the  equilit)rium  and  stability  of  toroidal 
plasmas  with  distributed  profiles.   Details  of  the  ideal  magneto- 
hydrodynamic model  are  given  by  Bernstein,  Frieman,  Kruskal,  and 

7 
Kulsrud.    The  total  energy 


E  =  fd\   {^PV^+Ib^  +j^] 


is  a  constant  of  the  equations  of  motion  assuming  that  fluid 
elements  expand  or  compress  adiabatically  during  any  motion  and 
assuming  (here)  that  there  are  no  electric  fields  or  body  forces. 
We  define  the  state  of  the  system  to  be  those  properties  (B  and  p) 
which  are  independent  of  the  velocity  v.   In  the  energy  variational 
principle,  we  wish  to  compare  the  potential  energy 


W 


^fd\   fiB^+^j  (26) 


of  only  those  states  which  are  accessible  to  each  other  by  a 
possible  motion  of  the  plasma.   The  condition  of  accessibility 
is  usually  built  in  to  the  variation  of  ¥  by  using  "constraints" — 
that  is  conditions  independent  of  velocity  which  are  to  be  the 
same  in  the  perturbed  and  unperturbed  states.   These  will  be 
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eniimerated  below.   It  has  been  shown  ^'   that,  subject 
to  these  constraints,  5W  =  0  is  equivalent  to  other  definitions 
of  static  equilibrium.   We  define  the  system  to  be  unstable  if 
and  only  if  any  other  accessible  state  has  a  minimized  potential 
energy  lower  than  that  of  the  equilibrium  being  considered. 

We  have  already  imposed  two  constraints  on  all  states  of  the 
system  to  be  considered  (to  apply  to  both  the  equilibrium  and 
perturbed  states): 

(1)  Nested  flux  surfaces  exist;  it  follows  that  B-VV  =  0: 

(2)  The  magnetic  energy  is  minimized  within  each  surface;  it 
follows  that  j.VV  =  0. 

In  order  to  minimize  the  rest  of  the  potential  energy  we 
now  demand  that 

(3)  The  thermodynamic  energy  is  minimized  within  each  surface.  3 
This  implies  that  the  plasma  is  to  be  in  thermodynamic  equilibrium 
within  each  surface  ;  the  pressure,  density,  temperature,  entropy, 
etc.  are  all  uniform  over  each  surface.   Hov/ever,  the  plasma  generally 
need  not  be  in  equilibrium  perpendicular  to  the  flux  surfaces.   In 
particular,  we  see  from  the  ideal  magnetohydrodynamic  equation  of 
motion 

p  dv/dt  =  -"^P  +  J  X  § 

that  constraints  (2)  and  (3)  Imply  that  there  is  a  component  of 
the  inertlal  term  p  dv/dt  only  perpendicular  to  each  flux  surface. 

After  satisfying  these  three  constraints,  the  potential  energy 
(26)  assumes  the  form 
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i  ^f    dV{|  4(V)Aij(V)  '^^(V)  +^]     .  (27) 

where  the  summation  convention  is  used  over  repeated  indicies 

i,j  =  1.2. 

Three  more  constraints  are  needed  to  assure  that  the 
perturbed  states  are  accessible  to  the  equilibrium  state.   These 
are : 

(4)  The  topology  of  each  field  line  within  the  perfectly- 
conducting  plasma  is  preserved  during  a  motion  of  the  fluid. 
Since  flux  is  carried  by  the  plasma,  it  follows  that  one  flux 
must  be  a  given  function  of  the  other  during  a  variation 

^2  =   ^2(^1)  (23) 

so  that  the  rotational  transform 

J:  (iy^)    =  27r  df^/di^-^ 

is  preserved  within  each  plasma  shell.   We  shall  take  the 
rotational  transform  as  part  of  the  data  which  describes  the 
equilibrium  state.   Eq.  (28)  implies  that  we  need  vary  only  one 
flux,  ^  =   'P-^{V)  . 

(5)  Entropy  is  conserved  during  the  motion  of  each  fluid  shell; 
that  is,  the  pressure  and  density  vary  according  to  the  adiabatic 
law. 

(6)  Fluid  mass  is  conserved  during  the  motion  of  each  fluid  shell. 
The  last  two  constraints  lead  to  a  relation  between  the  variation 
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of  pressure  and  the  variation  of  flux.   These  constraints  are  most 
easily  expressed  by  using  a  "convectlve  variation"  which  follows 
the  fluid  shells  rather  than  the  ordinary  variation  taken  at 
constant  volume. 

Consider  a  variation  of  the  flux  profile  as  shown  in  Fig.  (2). 
A  given  set  of  flux  surfaces  may  be  labeled  by  either  V  or  '//. 
Let  6  refer  to  the  variation  with  V  as  the  independent  parameter: 


5  -  5^  1^  (29) 


Let  A  refer  to  the  convective  variation  which  follows  a  given 
fluid  element,  where  ^  is  the  independent  parameter: 

^  =  ^V^  (30) 

Note  that  5V  =  0  and  A^  =  0.   We  can  graphically  determine  the 
inverse  function  relation  from  Fig.  (2): 


6^  =  -  fl  AV  (31) 


Then,  for  any  function  of  V  and  Tp   we  may  use  the  chain  rule 

_d_  ^  d£  a    S 
dV   dY'^       'Sy 

together  with  Eqs .  (29),  (30),  and  (31)  to  derive  the  identity 
relating  the  two  kinds  of  variation 
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A=5-^A  (32) 

The  shell  of  plasma  between  two  flux  surfaces  is  labeled 
by  the  fluxes  of  the  bounding  surfaces.   As  the  flux  profile 
is  varied,  the  plasma  moves  with  the  flux  so  that  each  plasma 
shell  retains  the  same  flux  labeling  but  may  change  its  volume. 
After  each  plasma  shell  comes  into  thermodynamic  equilibrium, 
we  may  express  the  adiabatic  compression  or  expansion  of  the 
plasma  by 


f=yf  (33) 


where  p  is  the  fluid  density  and  7  is  the  ratio  of  specific  heats 
(from  kinetic  theory,  7  =  l+2/(number  of  degrees  of  freedom)). 
In  order  to  calculate  Ap  we  need  to  use  the  conservation  of 
mass  within  each  fluid  shell,  constraint  (6), 

A(p  dV)  =0  .  (3^) 

Note  that  the  operators  A  and  d.  commute.   We  may  now  use  the 
last  three  equations  to  write 


^P=(i-r)^,|.A(rp|)  .  (35) 


Note  that  it  is  unnecessary  and  misleading  to  assume  that  the 
entropy  is  constant  over  the  plasma  profile.   We  have  assu.-ned 
only  that  the  entropy  and  mass  of  each  fluid  shell  remains 
constant  under  the  variation. 
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¥e  may  vary  the  flux  profile  and  the  inductance  either 
separately  or  both  simultaneously.   In  general,  we  must  also 
vary  the  total  current  profile  I^{V)    in  order  to  always  satisfy 
Eq.  (16).   However,  the  variation  of  1^{V)    never  appears  explicitly. 
The  variation  of  flux  and  inductance  separately  will  give  not  only 
all  the  possible  equilibriiwi  conditions  but  also  a  necessary 
and  sufficient  stability  condition.   That  Is^the  cross  term  in 
the  second  variation  with  both  flux  and  inductance  varied 
simultaneously  and  independently  does  not  affect  the  stability 
result. 

B.   Variation  of  Profile  with  Fixed  Inductance 

Conceptually  the  simplest  variation  of  the  minimum  potential 
energy  (27)  is  to  vary  the  flux  as  a  function  of  volume  while 
holding  the  inductance  as  a  function  of  volume  fixed.   That  is, 
we  rearrange  the  flux  labeling  of  the  surfaces  '/Without  changing 
their  shape.   We  take  the  flux  to  be  fixed  at  the  magnetic  axis, 
6^(0)  =  0,  and  at  the  perfectly  conducting  wall,  6^(C)  =  0. 
Depending  upon  which  flux  is  used,  one  condition  is  true  by 
definition  and  the  other  condition  may  be  relaxed  by  using  a 
free  boundary  variation. 

Using  the  fact  that  A-  •  is  a  symmetric  matrix,  we  have  for 

-'-J 

the  first  variation  of  potential  energy  (2) 


55=  +J    dV[^^^\.{^^P.)^+^v/{7-l)]    . 
Using  the  constraint  that  one  flux  is  a  function  of  the  other. 
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^^i  "  ^iV  ^^/^v  ' 


and  using  Eq.  (35)  to  relate  the  variation  of  pressure  to  variation 
of  flux,  and  integrating  "by  parts  with  fixed  end  conditons,  we 
find  that  y   drops  out  of  the  first  variation 


dj=  -  J      dV{  ^iv(Ai  j^j  v)  Y  +  dp/dV}  5^/^.^  . 


Our  first  equililDrium  equation  results  from  the  extremum  condition 
&^  =  0  for  arbitrary  flux  variations 

^iV^Aij  ^j^)v  +  dp/dV  =  0  (36) 


where 


^iV  ^  ^^i  (^)/<^^  • 


2 

The  second  variation  is  found  by  using  5^=0 


5^^.  =  y.  (&V^)^  -  ^*.  dTpbi/l'^ 


H^^/h  =-i^[7v^{Wh]/^+-iv^ky)/h^f/'/  ' 


After  some  algebra  we  find 


o 
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wnere 


6*  =   6^/-^   (assuming  ^(V)  7^  0) 

Is  an  arbitrary  test  function  with  &(1)(0)  =  6<i'(C)  =  0.   Note 
that  the  expression  in  square  brackets  is  positive  definite; 
hence  there  can  be  negative  contributions  from  only  the  last 
term  in  Eq.  (37).   The  final  minimization  for  this  part  of  the 
potential  energy  is  a  minimization  with  respect  to  the  test 
function  5<t)  which,  after  normalizing  54>,  leads  to  an  ordinary 
Sturm- Liouville  equation  for  54).   The  sign  of  the  lowest  eigen- 
value determines  the  stability  of  these  fixed  geometry  modes. 

In  the  special  case  of  straight,  coaxial  circular  cylinders 
(screw  pinch)  of  length  L  with  periodic  end  conditions  we  have 
the  following  inductance  matrix 

%1   ^    ■'^22    L  ^ 


A12  =-^21  =  °  • 


If  the  current  is  purely  longitudinal  and  concentrated  near 
the  axis,  it  is  possible  to  find  unstable  profile  distributions 
of  the  form 


^2(V)  ~  [aV/C-1]"-'-   for  a  »  1 


The  instability  is  a  breathing  mode  localized  near  the  axis. 
This  and  other  examples  will  be  discussed  in  future  papers. 
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C,   Incompressible  Variations 

For  the  second,  subclass  of  variations  we  imagine  that  the 
flux  surfaces  change  their  shape  but  that  each  plasma  shell 
preserves  its  volume  by  moving  to  a  new  surface.   Hence  we 
take  6^  =  0  and  therefore  Bp  =  0.   We  call  these  "incompressible 
variations"  because  the  fluid  shells  are  neither  compressed  or 
expanded.   However,  as  the  shells  change  shape,  individual  fluid 
elements  may  compress  or  expand  during  their  motion  within  the 
shell  as  they  reach  surface  equilibrium  (constraint  (3)).   Hence 
these  "incompressible  variations"  should  not  be  confused  with 
variations  of  an  incompressible  fluid.   The  latter  are  only  a 
subclass  of  the  former. 

We  label  the  flux  surfaces  by  their  volume  and  parameterize 

their  shapes  by  a  countable  number  of  functions  a  (V).   The 

functions  a  (V)  might  for  example,  be  the  ellipticity  of 

elliptical  surfaces,  or  the  relative  displacements  of  their 

centers,  or  bhe  amplitude  or  the  phase  of  a  helical  distortion, 

etc.   Since  the  specific  inductance  is  derived  from  first  order 

differential  equations  we  conclude  that  it  is  a  function  of  only 

the  zero  and  first  derivatives  of  a  (V) 

n^  ' 

^11  =Ai.(V;  a^(V),  a^'(V);  n  =  1,2,...) 


where  a  '(V)  =  d  a  (V)/dV;  no  derivatives  higher  than  first 
order  appear. 

We  vary  each  function  a  (V)  independently.   It  follows 
that  there  are  as  many  equilibrium  conditions  as  there  are  such 
functions;  each  has  the  form 


_Q^_ 
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^vAija/jV-  (^vAija'/jv)v  =  ^  (38) 

with 

^V  ^^ij/^a^'(V)    ^jv  Sa^(V)     |v=0   =  °    •  (39) 

Note  that  &a  =  0  at  V  =  C  because  the  wall  is  assumed  to  be 
n 

rigid  and  perfectly  conducting.  Since  the  total  volume  of  the 
system  remains  unchanged,  it  is  always  possible  to  take  BV'  =  0 
and  6p  =  0  at  V  =  0  and  at  V  =  C. 

The  equilibrium  conditions  given  by  Eq.  (33)  may  be  thought 
of  as  separate  extrema  for  each  one  parameter  family  of  surfaces 
being  considered.   If  the  class  of  surfaces  is  limited  —  for 
example,  to  elliptical  surfaces  with  variable  ellipticity  — 
Eq.  (33)  represents  a  constrained  equilibrium  —  as  opposed  to 
a  true  equilibrium  where  the  class  of  surfaces  is  limited  only 
by  smoothness  conditions.   The  concept  of  constrained  equilibria 
is  a  useful  tool  for  approximation.   But,  in  order  to  improve 
such  approximations,  we  need  to  characterize  the  class  of  sur- 
faces to  be  considered  by  a  complete  set  of  parameters.   The 
rest  of  this  work  is  independent  of  Eq.  (33);  all  other  results 
may  be  used  if  the  equilibrium  is  known  from  an  independent 
calculation. 

For  incompressible  variations  of  the  fluid  shells,  the 
stability  with  respect  to  infinitesimal  distortions  of  the  flux 
surfaces  can  be  determined  from  the  simple  form 

s"i  =  /^  dv  1  ^. ^  ^\.  i^.^  .  (40) 

It  has  been  verified  that  this  result  agrees  with  the  classical 
result  for  a  screw  pinch  where  an  arbitrary  infinitesimal  dis- 
tortion of  the  flux  surfaces  can  be  decomposed  into  a  spectrum 
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of  helical  distortions  with  k  ?^  0  or  m  ;^  0. 

We  may  consider  stability  with  respect  to  finite  distortions 
of  the  flux  surfaces  provided  p{V)    is  the  same  in  the  final  state 
as  in  the  equilibrium  state.   That  is,  we  replace  6^"  in  Eq.  (■^0) 
by  a  symbol  for  the  finite  difference.   We  are  therefore  able 
to  study  nonlinear  stability  with  respect  to  Incompressible 
variations.   We  are  not  able  to  find  a  simple  expression  for 
the  finite  variation  of  the  flux  profile  with  fixed  inductance 
because  the  conversion  between  convective  and  ordinary  variation 
relied  on  the  infinitesimal  nature  of  the  variations. 

D.   General  Variation 

In  general,  the  variations  5?//  and  &A  are  independent  of 
one  another  if  the  total  currents  1.  are  not  fixed  functions 
of  volume.   Eqs .  (36)  and  (38)  still  represent  all  the  equilibrium 
conditions  subject  to  our  constraints.   However,  cross  terms  are 
possible  in  the  second  variation  of  the  potential  energy 


C 


(41) 


1    ^■ 


+  ^v^^i,l(^v^*)v+^^v^^Ai,  V 


where 


5*  =  bf/f^,       f{V)    =  ^^(V) 


If  either  the  fixed  inductance  variation  or  the  incompressible 
variation  were  unstable  we  would  need  to  go  no  further  than  to 
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consider  the  special  second  variations  given  by  Eqs .  (37) 
or  CJO).   However,  If  both  special  variations  were 
positive,  we  would  need  to  examine  the  effect  of  the  cross 
term  In  Eq.  Cll)  In  order  to  find  a  sufficient  as  well 
as  necessary  condition  for  stability. 


E.   Equivalence  with  the  Classical  Energy  Principle 

We  shall  show  that  the  energy  principle  that  has 
been  derived  here,  using  constraints  (1)  through  (6),  is 

equivalent  to  the  partially  minimized  ordinary  energy 

7 
principle  when  used  for  the  purpose  of  distinguishing 

stability  from  instability.   However,  growth  or  decay  rates 
cannot  be  calculated  from  our  variational  principle  di- 
rectly. 

7 
Bernstein,  Frieman  ,  Kruskal  and  Kulsrud  point  out 

in  section  4b  that  it  is  sufficient  to  normalize  only  that 

component  of  the  displacement  S which  is  orthogonal  to 

the  p   =  constant 
o 


-26- 


surfaces  in  order  to  distinguish  stability  from  instability. 
The  condition  &W  be  minimum  yields  an  Euler  equation  with  three 
components.   With  only  the  orthogonal  component  of  |  normalized, 
the  two  components  of  the  Euler  equation  in  the  surface  are,  to 
first  order  in  |, 

j*Vp  =  0   and  B.Vp  =  0  . 

If  we  assume  the  existence  of  well  defined  flux  surfaces  (constraint 
(1))  it  follows  from  these  equations  that  p  is  constant  on  flux 
surfaces  (constraint  (3))  and  that  J»VV  =  0  (constraint  (2)). 
Hence,  for  infinitesimal  perturbations  we  have  equivalence  of  the 
energy  principles  in  the  sense  that 

min  &W   ^  0  ±5-  6  5  ^  0  . 


IV   Special  Calculations  of  Specific  Inductance 

It  seems  to  be  possible  to  find  exact  formulas,  in  terms 
of  one  integration,  for  the  specific  inductance  of  any  nested 
set  of  smooth  flux  surfaces  with  an  ignorable  coordinate.   We 
shall  derive  such  formulas  for  a  straight  cylinder,  an  axisy- 
mmetric  torus,  and  a  straight  helically  symmetric  system,  each 
with  an  arbitrary  cross  section.   For  the  first  two  cases  we 
need  to  use  only  Ampere's  law,  without  appealing  to  the  general 
methods  derived  in  section  II.   For  the  case  of  helical  symmetry 
we  use   a  method  suggested  by  Grad.   That  is,  we  find  two 
simple  solutions  for  the  two  potentials  of  Eq.  (l8)  and  then 
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determine  their  periods. 

In  the  case  of  a  straight  cylinder  and  an  axisymmetric 
torus,  it  is  possible  to  perform  the  integration  analytically 
when  the  cross  sections  are  circles  with  relative  displacements 
"between  their  centers.   By  collecting  these  exact  solutions  we 
may  begin  to  develop  intuition  about  the  specific  inductance  in 
general. 

A.   Straight  Cylinder 

We  consider  straight  cylindrical  flux  surfaces  with  arbitrary, 
smooth  cross  sections  given  by  V  =  V(x,y).   If  the  only  component 
of  current  is  parallel  to  the  ignorable  z  coordinate  then  the 
vector  potential  has  only  a  z  component  B  =  V  x  A  =  VA(x,y)  x  z 
and  all  fl\ix  is  aximuthal.   If  we  set  A  =  0  along  a  reference 
line  parallel  to  the  z  axis,  we  can  use  Stoke 's  theorem  on  the 
definition  of  the  flux  [Eq.  (l) ]  to  show 

Y(x,y)  =  A(x,y)  L 

where  L  is  the  periodic  length  of  the  cylinder.   This  flux  is 
a  function  of  V  alone 

¥(x,y)  =  ^-L(V(x,y))  . 

These  equations  are  used  in  Ampere's  law  Eq.  (3)   to  find  the 
relation  between  longitudinal  current  and  azimuthal  flux 
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J,2   =   f    di.B  =  6      di.i  VY  X  z 


'-'P        L-p 


=  ^  qT  d^'VV(x,y)  X  z  ^^(V) 


"1 

Hence,  the  azimuthal  specific  inductance  is 

A22(V)  =i/  ^|-VV(x,y)  X  z  (42) 

^2 

where  Cp  is  any  closed  contour  on  the  flux  surface  wrapping  once 
around  the  cylinder. 

For  the  calculation  of  the  longitudinal  specific  inductance, 
A-,-,  we  place  contour  C-,  along  the  ignorable  z  coordinate  with 
the  return  element  along  a  reference  line.   From  Ampere's  law 
we  find  a  relation  bet'.jeen  the  azimuthal  current  and  the 
longitudinal  field 


T   =  B.     L 
-J-.az    long 

Since  J   is  independent  of  the  choice  of  C-,  contour  taken  on 

any  given  flux  surface,  we  conclude  that  B-,     is  uniform  over 
•'to  '  long 

each  flux  surface,  regardless  of  the  cross  sectional  shape.   If 
do-,  is  the  area  of  the  cut  surface  S-,  between  differentially 
close  flux  surfaces,  we  use  dV  =  L  dS^  to  show 

d^p         dS,  P 

-^   =  B,    -:^,   =  B,    /L  =  T   A 
dV    long  dV    long     -^  az 


Hence,  the  longitudinal  specific  inductance 
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,2 


All  =  ^  ("^3) 


is  independent  of  the  cross  sectional  shape  of  the  flux  surfaces. 
Note  that  there  is  no  mutual  inductance  for  straight  cylinders. 

One  of  the  simplest  non-trivial  applications  of  Eq.  (42)  is  to 
the  case  of  non-concentric  circular  cylinders.   If  r  is  the 
radius  of  such  a  surface  and  a(r)  is  the  displacement  of  its 
center  along  the  x  axis,  then  these  surfaces  may  be  represented 

T3y 


V/ttL  =  r^  =  y^+(x-a(r)) 


This  form  may  be  used  in  Eq.  (1)  to  find 

Vl-   (a'(r))^ 

If  a(r)  =  0,  we  have  the  specific  azimuthal  inductance  for 
concentric  circular  cylinders 


A22=T^  (45) 


This  may  be  derived  from  the  ordinary  vacuum  field  inductance, 
which  is  proportional  to  ln(r),  by  first  applying  d/dV  and  then 
taking  the  reciprocal.   Returning  to  the  nonconcentric  result 
(3),  we  note  thatApp  ~*  oo  as  |  da/dr  |  -*  1.   Geometrical  analysis 
indicates  that  differentially  close  flux  surfaces  touch  and  cross 
each  other  when  | da/dr |  >  1.   Intuitively,  then,  the  specific 
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inductance  behaves  like  the  capacitance  between  close  surfaces 
and  it  becomes  infinite  as  they  touch. 

B.   Axisymmetric  Torus 

If  we  use  cylindrical  coordinates  (R,  0,  z)  and  take  <t>   as 
the  ignorable  coordinate,  a  calculation  similar  to  the  above 
yields  the  specific  inductance  for  axisymmetric  toroidal  flux 
surfaces  with  arbitrary  cross  sections  and  relative  displacements 

A22  "  /  '^^'^^  ^   V(2-/rR)  (46) 

^2 


All  =  2-dV//     f-  .  (^7) 

^1 

Again  the  matrix  is  diagonal. 

As  a  special  case,  consider  symmetric  toruses  with  circular 
cross  sections   of  radius  r  with  their  centers  displaced  in  the 
plane  of  the  torus  by  a(r)  relative  to  the  magnetic  axis  at 
R  =  R  ,-  positive  values  of  a(r)  represent  a  displacement  inward 
from  the  magnetic  axis  toward  the  center  line.   Such  toruses  may 
be  represented  by  the  equation 


V  =  2Tr^(R  -a)r^  =  27r^(R  -a)  [  z^  +  (R-R  +a)^]  . 


If  we  let  D  =  da/dr  and  e  =   r/R  (the  aspect  ratio),  the  exact 
results  are 


A    -    2V/R      f 1_  + D/§ 1  (l^Q) 

/l-E^   ^(l+£D/2)2-D^ 


22    l+D/e  +  eD72 
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s/i:7 

2     2 
where  (l+eD/2)   >  D  is  the  condition  that  no  surfaces  touch 

each  other.   Note  that  the  inductances  are  finite  at  l+D/e+eD/2  =  0, 

g.   Helical  Symmetry 

For  systems  topologically  equivalent  to  a  periodic  cylinder, 
the  computations  are  most  easily  carried  out  if  we  retain  the 

/\   /\   >^ 

unit  vectors  of  cylindrical  coordinates  (z,9,r)  but  change 
variables  to  (z,0, R)  where  R  =  R(z,  9,  V).   The  z  and  6   components 
of  Eq.  (l8)  are 


-^^^e  =   (i^W)*z-^o^z*0 


^  Q  V  =   -RflR  *  +(1+R  ^)*fl 
Ry.  ^  z     0  z  z   ^    z    '    Q 


(50) 


where 


Q  =  1  +R^^  +Rg^/R^  . 


The  r  component  of  Eq.  (l8)  is  a  linear  combination  of  these 
equations.   With  these  conventions,  the  periods  are 


6.9   <t)^ 
O       ~     "    ^  o 


ii  =  /,  ^^  \'  I2  =  X 


27r        0     r^ 

d0  Vq,   -//^  =  /   , 
o  "^  o 


(51) 


Any  surface  with  helical  symmetry  may  be  represented  by 
R  -   R(^,V)  where  ^  =  mS  +kz.   If  the  surfaces  were  straight 
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circular  cylinders,  where  V  =  ttR^L,  the  solutions  to  Eq.  (50) 
would  oe 


V  =  C^z  +(l^e 


<J)   =   +2\AC^e  -  GTTLCgZ    . 


A   solution  is   found   for  the  helical   system  by  adding   v(^,V) 
and   *(C,V")    to   these   solutions.      From  Eq.    (50),    these  new 
functions   satisfy  the  ordinary  differential  equations 


(m^+:^-R^)   -||  =  k[2VRR^-R^]C^ +m[27rLRR^-l]C2 


(m^  +A^)   ^  =  m[2V-  RQ/R^]C-^  +k[RQ/R^  -  27rLR^  ]C2 


The  results  of  inserting  the  inegrals  of  these  equations  [Eq.  (51)] 
into  Eq.  (l6)  and  solving  for  A^  .•  are 

f^        =  L^ { 1  +  PP  +  a ( M+N)  +  am^  ( M+N")  2/(  i  -  aN  -  nrM)  ] 

y\_22   =  "X^tl  +n'P  +aN  +a6N^/(l  -  aN  -  m^M) }  (52) 


A  n  2  "  -^Pl   "  2Vmk(  P  -  aN( N+M) /(I  -  aN  -  m^M) } 


Here  we  have  represented  the  surfaces  by 


ttR^L  =  (1  +a(^,V))V 
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and  we  have  used  the  notations 


Va,, 

M  .  -  <-^>,   N  .  -  < ^> 

a+Pa  a+Pa 

T  a.^  V^.JV8i)„ 


l+a 

a+Pa 

<.. 

..> 

_  r^^d^ 

a 

=  m^+P 

?   = 

:  Vk^/frL   . 

No  approximations  have  been  used   to  derive  these  results  from 
Eqs.  (50),  (51),  and  (l6). 
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V.   Conclusions 

An  essential  step  In  this  formulation  of  the  energy 
principle  Is  the  minimization  of  potential  energy  before 
comparing  states  with  a  variational  analysis.   It  Is  shown 
that  a  unique  magnetic  field  minimizes  the  magnetic  energy 
for  a  given  set  of  nested  flux  surfaces  and  flux  profiles. 
For  this  minimizing  field,  the  specific  magnetic  energy  can 
be  factored  Into  geometric  and  profile  dependent  matrlcles 
as  In  Eq.  (17).   The  geometry-dependent  matrix,  called  the 
specific  Inductance  A-.>  Is  defined  by  Eq.  (l6).   It  can 
be  represented  by  Eq.  (22)  or  Eq.  (25)  If  and  only  If  we 
use  the  coordinate  system  defined  by  Eqs .  (21)  and  (23). 
Alternatively,  If  any  basis  pair  of  total  currents  and 
differential  fluxes  Is  known  for  a  given  set  of  flux  sur- 
faces, A...  may  be  computed  directly  from  Eq.  (I6). 

The  equilibrium  conditions  given  by  Eqs.  (36)  and  (38) 
result  from  extremlzlng  the  potential  energy  with  respect 
to  flux  profile  and  surface  shapes.   It  Is  necessary  and 
sufficient  for  stability  that  the  second  variation  given 
by  Eq.  (41)  be  positive.   Eq.  (37)  represents  the  special 
case  of  a  variation  which  leaves  the  surface  shapes  un- 
changed (fixed  Inductance).   Eq .  (40)  represents  a  varia- 
tion which  leaves  the  flux  profile  unchanged  (Incompressible 
variation).   The  latter  can  be  easily  generalized  to 
consider  finite  distortion  of  the  flux  surface  shapes  and 
hence  nonlinear  stability.   It  Is  shown  In  section  III-E 

that  the  above  stability  criterion  Is  equivalent  to  that 

7 
of  the  classical  energy  principle. 
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There  are  a  number  of  promising  applications  for  this 
new  formulation.   The  one  dimensional  character  of  the 
integrations  and  the  simple  decomposition  of  the  integrands 
makes  this  formulation  especially  suitable  for  computer 
studies  of  complicated  geometries.   The  use  of  restricted 
families  of  surface  shapes  (circular,  elliptical,  etc.)  for 
which  there  are  simple  calculations  of  the  specific  inductance 
is  an  aide  to  the  analytic  approximation  of  equilibria  and  the 
formulation  of  necessary  conditions  for  stability.   The  coordi- 
nate system  defined  by  Eqs .  (21)  and  (23)  provides  an  alterna- 
tive to  Hamada  coordinates  for  the  representation  of  stability 
criteria  using  a  general  perturbation  analysis.   Finally,  the 
study  of  finite  incompressible  variations  may  be  a  useful 
approach  to  the  study  of  nonlinear  stability. 
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Appendix 

The  following  two  derivations  may  be  found  in  Grad  and 
Rubin^  and  in  Kruskal  and  Kulsrud.    They  are  presented  here 
without  any  reference  to  equilibrium  and  with  a  redefinition 
of  the  coordinate  conventions. 

I.    The  Periods  of  v 

Given  nested  toroidal  surfaces  V  =  V{9      9     r)    and  the 
conventions  illustrated  in  Fig.  1  we  employ  Stokes'  theorem 
to  show 

cf   di.v(^)v  =  r      di-y^^Ke   -  r      diM^Ke 
=  -/  d^.v(2)v02  =  f   ds -|^^.vx(v(2)ve2). 

The  integrand  is  transformed  with  vector  identities,  including 

VV.V  X  V^^^v  =  0 
to  show 


VV.V  X  (V^^^vSg)  =  '^'{9^7V  X  v(2)v) 


¥e  then  use 

■^  S^'JWT  '"   ~  dvjy  ax... 

and  the  divergence  theorem  over  the  following  differential  volume; 
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>  dS^ 


only  the  surface  at  9     =   1  survives;  hence 


/  di'V(^)v  =  -§^  f    dS-VV  X  v(^)v 


The  same  steps  may  be  followed  to  show  that 


(C    d^.v(^)v  =  -  A  r  dS.VV  X  v(^)v  . 


These  equations  imply  Eq.  (9) • 


II   Application  of  Hodges '  Theorem 


Given 


B  = 


=  VV  X  V  [0j_^2"  °2^1  "^^^ 


VV.V  X  B  =  0 


A  single  valued  on  V  =  V{9^,9^,r) 


Show  that 
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fd\      B^  =Xr   r  dV  i.(V)  ^  (V) 


starting  with 


jd^x  B^  =  ^d^x   B.VV  X  V[0^^2  -^2^1  "^^^ 


w 


e  use  vector  Identities,  Eq.  (l4)  and  Stokes'  theorem  to  show 


jd^x  B^  =ydv  ^ds  1^.  V  X  {[eJj^-e,J^  +  },]B) 
o 

=  /dV  cT  di.B[e^^2 -^2^1  +^] 


^"l 


With  the  above  illustrated  contours,  we  show 


"^    2 


•      / 


'd-^    B^  =  AvLV'g  5^   di.B  +  ^^  c^   d^.B] 


and  hence  the  result. 
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Figure  Captions 

Pig.  1    The  orientation  of  the  contours  C,  and  Cp  and  cut 
surfaces  S-,  and  Sp  are  Illustrated  for  a  given  toroidal  flux 
surface  V  =  V{e^,9^,r). 

Fig.  2    The  relation  between  the  ordinary  variation  of  flux 
as  a  function  of  volume  5^  and  the  convective  variation  of 
volume  as  a  function  of  flux  AV  is  illustrated. 
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